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The inclusive mesons production in proton-nucleus collisions in the near threshold and 
subthreshold energy regimes is analyzed with respect to the one-step {pN —> K^YN, Y = A, S) and 
the two-step {pN NNtt, NN2tt-, ttN —> K'^Y) incoherent production processes on the basis of an 
appropriate new folding model, which takes properly into account the struck target nucleon removal 
energy and momentum distribution (nucleon spectral function), extracted from recent quasielastic 
electron scattering experiments and from many-body calculations with realistic models of the NN 
interaction. Comparison of the model calculations of the K'^ double differential cross sections for 
the reaction p + C^^ at 1.2, 1.5 and 2.5 GeV beam energies with the existing experimental data from 
the SATURNE experiment is given, illustrating both the relative role of the primary and secondary 
production channels at considered incident energies and those features of the cross sections which 
are sensitive to the high momentum and high removal energy part of the nucleon spectral function 
that is governed by nucleon-nucleon short-range and tensor correlations. In-medium modifications 
of the available for pion and hyperon production invariant energies squared due to the respective 
optical potentials are needed to account for the experimental data on production in the energy 
region under consideration. 


INTRODUCTION 

Extensive investigations of the production of K'^ mesons in proton-nucleus reactions [1-9] at incident energies 
lower than the free nucleon-nucleon threshold have been carried out in the last years. The study of inclusive and 
exclusive subthreshold K'^ production in pA interactions is planned in the near future at the accelerators COSY- 
Jiilich [10] and CELSIUS [11]. Because of the high K~^ production threshold (1.58 GeV) in the nucleon-nucleon 
collision and the rather weak K'^ rescattering in the surrounding medium compared to the pions, etas, antiprotons 
and antikaons, from these studies one hopes to extract information about both the intrinsic properties of target nuclei 
(such as Fermi motion, high momentum components of the nuclear wave function, clusters of nucleons or quarks) 
and reaction mechanism, in-medium properties of hadrons. However, in spite of large efforts, subthreshold kaon 
production in proton-nucleus reactions is still far from being fully understood. Thus, measured total production 
cross sections [1] in the range of proton energies of 800-1000 MeV were described in the framework of the respective 
folding models based on both the direct mechanism [2,6] of K'^ production [pN —> K^AN) and on the two-step 
mechanism [3-5] associated with the production of kaons by intermediate pions {pNi irNN, ttN 2 A) using 

different parametrizations for internal nucleon momentum distribution [6, 12, 13]. The same experimental data could 
be well explained also in the modified phase space model [7]. In the above folding models only the internal nucleon 
momentum distribution has been used and the off-shell propagation of the struck target nucleon has been neglected 
or has been taken into account the most crudely, but it could be significant in processes that are limited by phase 
space such as the threshold heavy mesons production. As is well known [13-21], the off-shell behaviour of a bound 
nucleon is described by the nucleon spectral function P(pt, E), which represents the probability to find in the nucleus 
a nucleon with momentum pt and removal (binding) energy E and contains all the information on the structure of a 
target nucleus. The knowledge of the spectral function P{pt,E) is needed for calculations of cross sections of various 
kinds of nuclear reactions. In particular, it has been widely used earlier for the analysis of inclusive and exclusive 
quasielastic electron scattering by nuclei [13-26]. It was found that even at very high momentum and energy transfer 
in this scattering, the target nucleus cannot be simply described as a collection of A on-shell nucleons subject only to 
Fermi motion, but the full nucleon momentum and binding energy distribution has to be considered. Only recently [8, 
9] the binding energy of the struck target nucleon has been properly taken into account in calculating the subthreshold 
kaon production in pA collisions. Debowski et al. [8] analyzed their own double differential cross sections data for 
the K~^ production in p C reactions at 1.2, 1.5 and 2.5 GeV beam energies in the framework of the first chance 
collision model with the nucleon spectral function taken from [20, 27]. It was shown that at subthreshold incident 
energies the production cross sections are underestimated significantly by calculations assuming only first chance 
collisions. The same nucleon spectral function has been employed by Sibirtsev et al. [9] for the description of the 
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measured [1] total cross sections for production from p + C collisions in the framework of the two-step model. 
The claims have been advanced in [9] that within the spectral function approach the two-step production mechanism 
with an intermediate pion also dominates at subthreshold energies as in the folding models [3-5], but the calculated 
total cross sections for production are underestimated below about 920 MeV in contrast to the folding model 
where a better description of the experimental data is achieved as well as the contribution to the calculated total cross 
sections from the high momentum and high removal energy part (correlated part) of the nucleon spectral function 
is small compared to that calculated with the single-particle (uncorrelated) part of the one. However, in order to 
fully clarify the role played by nucleon-nucleon correlations induced by realistic interactions in the subthreshold kaon 
production in pA collisions as well as to get a deeper insight into the relative role of the primary and secondary 
reaction channels, it is necessary to properly analyze within the same approach, based on nucleon spectral function, 
the existing experimental data both on total [1] and differential [8] kaon production cross sections. 

In this paper we have performed such description of the total and differential K'^ production cross sections from 
pC reactions in the near threshold and subthreshold energy regimes using the appropriate new folding model for 
primary and secondary production processes, which takes into account the struck target nucleon removal energy and 
momentum distribution. 


I. DIRECT K+ PRODUCTION PROCESS 


Apart from participation in the elastic scattering an incident proton can produce a directly in the first inelastic 
pN collision due to nucleon Fermi motion. Since we are interested in a few GeV region (up to 2.5 GeU), we have 
taken into account [28] the following elementary processes which have the lowest free production thresholds: 

p + N ^ K+ + A + N, (1) 


p + N ^ K+ + T, + N. 


( 2 ) 


Neglecting the kaon rescatterings in the nuclear medium [29], we can represent the invariant inclusive cross section of 
AT'*" production on nuclei by the initial proton with momentum po as follows [13, 23, 29]: 


Ek+ 


d^K+ 


= /v[A]x 


( 3 ) 


where 


Ek+ 


dcFpN^K+KN{VOiVK+) 

dpK+ 


+ ( Eii+ 


dgpjv^jc+sjv(Poi Pj^tx) 
dVK+ 



dc^pN —> 


ic+yAf(Poi Pic+) 
dpK+ 


j J P{pt,E)dptdEx 


( 4 ) 


Ek+ 


dapN^K+ yn{Vs, Pk+ ) 
dpK+ 


Here, Ex+dapf^^x+YNiV^i PK+)/dpK+ are the free invariant inclusive cross sections for production in reactions 
(l)-(2); P{pt,E) is the nucleon spectral function normalized to unity; pt and E are the internal momentum and 
removal energy of the struck target nucleon just before the collision; Pk+ and Ej^+ are the momentum and total 
energy of a K'^ meson, respectively; /v[A] is the effective number of nucleons for the pN K'^YN reaction on nuclei 
1^; s is the pN center-of-mass energy squared. The expression for s is: 

s = (A'O + Et)'^ — (po + Pt)^, (5) 


'^Calculations show that in the case of Gaussian nuclear density for target nucleus the effective number of nucleons is 
equal to 6.9. 
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where Eq and Et are the projectile’s total energy, given by Eq = a/pq + rnj^ (niN is the rest mass of a nucleon), 
and the struck target nucleon total energy, respectively. Taking into account the recoil and excitation energies of the 
residual (A — 1) system, one has [23]: 

Et = MA- ^/{-pt)^ + {MA-mN + E)^, ( 6 ) 

where Ma is the mass of the initial target nucleus. It is easily seen that in this case the struck target nucleon is 
off-shell. In Eq. (3) it is assumed that the meson production cross sections in pp- and pn-interactions are the 
same [8,28] as well as it is disregarded any difference between the proton and the neutron spectral functions [18]. In 
our approach the invariant inclusive cross sections for production in the elementary processes (l)-(2) have been 
described by the three-body phase space calculations normalized to the corresponding total cross sections [28]. 

The nucleon spectral function, P{pt, E), is a crucial point in the evaluation of the subthreshold production of any 
particles on a nuclear target. When ground-state NN correlations, which are generated by the short range and tensor 
parts of realistic NN interaction, are considered, the spectral function P{pt,E) can be represented in the following 
form [16-21]: 

Pipt,E)=Po{pt,E)+Piipt,E), (7) 

where Pq includes the ground and one-hole states of the residual (A — 1) nucleon system and Pi more complex 
configurations (mainly lp-2h states) that are arise from 2p-2h excited states generated in the ground state of the 
target nucleus by NN correlations. In the absence of NN correlations (within the Hartree-Fock approximation) 
function Pi{pt,E) = 0 and the Hartree-Fock spectral function [16-21] is recovered. Thus, the NN correlations 
deplete (20 ± 5%, [16, 17, 21, 26, 30]) states below the Fermi sea and populate states outside the Fermi sea. For 
K'^ production calculations in the case of target nucleus reported here we have employed for the single-particle 
part Po{pt,E) of the nucleon spectral function the following relation [16-18]: Po{pt,E) = 0.8P^^^'>{pt, E) with 
{pt, E) being the harmonic oscillator spectral function (see, e.g., [24]) in which the s- and p-shell nucleon 
momentum distributions were taken from [29] and binding energies of 34 and 16 MeV [24] for the s and p shells, 
respectively, were used. 

Results of the many-body calculations with realistic models of the NN interaction in few-body systems, complex 
nuclei and nuclear matter [16, 17, 19-23] as well as experimental data for nucleon momentum distribution n(pt) = 
J P{pt, E)dE obtained by the y-scaling analysis [18] show that: 

i) the momentum distribution n(pt) for pt > 2fm~^ is a several orders of magnitude larger than the predictions 
from mean-field calculations and, moreover, the high momentum behaviour of n(pt) (at pt > 2fm~^) is similar 
for nuclei with mass number A > 2; 

ii) the behaviour of the nucleon spectral function at high values of pt and E is almost entirely governed by Pi (pt, E) 
as well as there is a strong coupling between the high momentum components and the high values of the removal 
energy [E ~ pl/2mN) of a nucleon embedded in the nuclear medium. 

Starting from such observation and analyzing the perturbative expansion of the NN interaction and the momentum 
distribution for potentials decreasing at large values of pt as powers of pt, Frankfurt and Strikman have argued 
[31] that the structure of the nucleon spectral function at high values of nucleon momentum and removal energy 
should be generated mainly by the ground-state two-nucleon configurations with large relative (> 1.5 but 

low (<1.5 fm~^) center- 0 f-mass momenta. Basing on this assumption, the convolution model for the correlated part 
Pi{pt,E) of the nucleon spectral function in the case of any value of the mass number A has been proposed [19, 21] 
according to which function Pi{pt, E) is expressed as a convolution integral of the momentum distributions describing 
the relative and center-of-mass motion of a correlated NN pair embedded in the nuclear medium. It was shown [19, 
21] that this model satisfactorily reproduces the existing spectral functions oi He^, He'^ and nuclear matter obtained 
by means of many-body calculations with realistic models of the NN interaction. An inspection of the convolution 
formula (53) from [21] for the spectral function Pi{pt,E) leads to the simple analytical expression for the Pi{pt,E) 
proposed in [17] (formula (7)). This expression for the quantity Pi{pt, E) (in which the correlated part ni(pt) of the 
n(pt) in the case of target nucleus was taken from [16], the peak position and the FWHM were taken from [21]) 
was used in our calculations of production in pC collisions. 

Let us consider now the two-step K'^ production mechanism. 


II. TWO-STEP K+ PRODUCTION PROCESS 

Kinematical considerations show that in the bombarding energy range of our interest (< 2.5 GeV) the following 
two-step production process may not only contribute to the production in pA interactions but even dominate 
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[1, 3-5] at subthreshold energies (< 1.5 GeV). An incident proton can produce in the first inelastic collision with an 
intranuclear nucleon also a pion through the elementary reactions: 

p + + N + (8) 


p + Ni^ N + N + 2 tt. (9) 

We remind that the free threshold energies for these reactions respectively are 0.29 and 0.60 GeV. Then the interme¬ 
diate pion, which is assumed to be on-shell, produces the kaon on a nucleon of the target nucleus via the elementary 
subprocesses with the lowest free production thresholds (respectively, 0.76 and 0.89 GeV): 

tt + N2^K+ +K, ( 10 ) 


■k + N2^ K+ + T., 


( 11 ) 


provided that these subprocesses are energetically possible. In order to calculate the K'^ production cross section for 
pA reactions from the secondary pion induced reaction channels (10), (11) we fold the momentum-energy-averaged 
inclusive differential cross section for pion production in the reactions (8), (9) (denoted by < dapN^T^xiVo^ P 7 r)/dp 7 i- >) 
with the momentum-energy-averaged inclusive invariant differential cross section for AT'*' production in these channels 
(denoted by < Ex+daTrAr^K+xiP-K, PK+)/dpK+ >) and the effective number of NN pairs per unit of square (denoted 
lia /.[A, J-M, (p.), ^.]) i.; (see! also, [29]): 


Ek+- 


da 


(sec) 

pA^K+X 




dpK^ 


( 12 ) 


d<TpX — >7tX (PO: Ptt) 
dp^ 


Ek+ 


d<J-KN^K+x{P'A^PK+) 

dpK+ 


Here, p,r is the momentum of a pion. The quantity (p,r),'d,r]) in Eq. (12) is determined by the 

formula (39) from [29] in which one has to make the following substitutions: NZ ^ A^, rj ^ TT, and ^ ai%. It 
should be noted that the symbol f dp,r in Eq. (12) stands for summation over the pion isospin and integration over 
the pion momentum space. 

In our method the differential cross sections for pion production in the elementary processes (8), (9) have been 
described by the three- and four-body phase space calculations normalized to the respective total cross sections [32]. 
The invariant inclusive cross sections for K~^ production in the elementary reactions (10), (11) have been described 
by the two-body phase space calculations [29] normalized to the corresponding total cross sections [33] 

and ct-kn^k+t. [34]. We have taken into account in the calculation of the Ar+ production cross section (12) from the 
two-step processes (8)-(ll) also the following medium effects the influence of the average momentum-dependent 
optical potential on the incoming proton (denoted by V^), on produced in the reactions (8), (9) nucleons (denoted by 
Vn), (10), (11) hyperons (denoted by Vy) □ . The available in this case for the pion and hyperon production invariant 
energies squared of the pN- (s) and 7rA^-systems (si) can be written as follows (cf. with Eq. (5)): 


^ — (Eq ~ 2bjv + Et)^ — (pq -|- pt)^. 


(13) 


Si — (Tf-TT ~ Vy -\- Et)"^ — (p,r + Pt)^) (14) 

where Eq, pg are the in-medium total energy and momentum of the initial proton (Eg = Eq — Vq, p^ = Pq — 2mivVb — 
Vq); Ett is the total energy of the intermediate pion. 

Now, let us discuss the results of our calculations in the framework of approach outlined above. 


^We have neglected the analogous medium effects in calculating the E+ production cross section (3) in the one-step processes 
(1), (2), since the contributions from them to the invariant energy squared s compensate each other. 

®The pion and kaon potentials in a nuclear medium were set to zero [35]. 
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III. RESULTS AND DISCUSSION 


The comparison of the results of our calculations by (3)-(7), (12)-(14) with the experimental data [8] for the double 
differential cross sections for the production of A+ mesons at an angle of 40° in the interaction of protons with 
energies of 1.2, 1.5 and 2.5 GeV with nuclei is given in Figs.1-3. It is seen that: 

1) the two-step production mechanism clearly dominates at subthreshold beam energies in accordance with 
earlier calculations [3-5] based on the folding model with one-dimensional internal nucleon momentum distribu¬ 
tions and it is of minor importance at 2.5 GeV proton beam energy at which mesons are produced almost 
entirely in first chance collisions (1), (2) in line with our earlier calculations [28]; 

2) the main contribution to the K'^ production in the two-step reaction channels (8)-(ll) both at subthreshold and 
above the free NN threshold incident energies considered here comes from the use only of the uncorrelated part 
Poipt, E) of the nucleon spectral function in the calculation of the corresponding momentum-energy-averaged 
differential cross sections for pion and kaon production; 

3) the K~^ production cross section in the two-step reaction channels (8)-(ll) coming from the use of the uncor¬ 
related and correlated parts of the nucleon spectral function in the calculation of momentum-energy-averaged 
differential cross sections for pion and kaon production, respectively, is approximately the same as that obtained 
from the use of the correlated and uncorrelated parts of one in the analogous calculation; 

4) the contribution to the production in the same reaction channels only from the correlated part Pi{pt,E) 
of the nucleon spectral function is an order of magnitude smaller than the contributions from the use of the 
uncorrelated (correlated) and correlated (uncorrelated) parts of one in calculating the momentum-energy- 
averaged differential cross sections for pion and kaon production, respectively; 

5) the contribution to the production from the secondary reaction channels (10) and (11) with a A and a 
S particles in the hnal states are comparable at 2.5 GeV beam energy, whereas at subthreshold energies the 
secondary production process (10) is more important than the (11) one; 

6) our full calculations (the sum of results obtained both for the one-step (1), (2) and two-step (8)-(ll) reaction 
channels, solid lines) reproduce the experimental data only if we take into account in the two-step production 
processes (8)-(ll) the in-medium modifications according to Eqs.(13), (14) of the available for pion and hyperon 
production invariant energies squared of the pN- and 7rA-systems. 

Finally, it should be noted that our model calculations for the two-step production processes (8)-(ll) with the use 
only of uncorrelated part of the nucleon spectral function as well as the same in-medium modifications of the available 
for pion and hyperon production invariant energies squared as in above reproduce also the measured [1] total 
production cross sections in pG collisions at proton energies above about 900 MeV but underestimate the data at 
lower bombarding energies, what indicates the need for other reaction channels or multinucleon correlations to explain 
these data close to the absolute production threshold. 

Taking into account the considered above, we conclude that the high momentum and high removal energy part of 
the nucleon spectral function which is generated by ground-state two-nucleon short-range correlations apparently 
cannot be studied in the inclusive subthreshold kaon production in pA reactions. As it was shown in [29], the 
7r+ + A —>■ K~^ + X reaction in the subthreshold regime seems to be quite promising for this aim. 


IV. SUMMARY 

In this paper we have calculated the total and differential cross sections for K'^ production fromp-|-C^^ reactions in 
the near threshold and subthreshold energy regimes by considering incoherent primary proton-nucleon and secondary 
pion-nucleon production processes within the framework of an appropriate folding model, which takes properly into 
account the struck target nucleon momentum and removal energy distribution. The comparison of the results of our 
calculations with the existing experimental data [1, 8] was made. It was found that the in-medium modifications of 
the available for pion and hyperon production invariant energies squared due to the respective optical potentials are 
needed to account for the experimental data both from Koptev et al. [1] and from the SATURNE/GSI collaboration 
[8]. It was shown also that the two-step production mechanism clearly dominates at subthreshold energies in line 
with earlier calculations [3-5] based on the use of the one-dimensional momentum distributions of nucleons in the 
target nucleus. Our investigations indicate that the main contribution to production in the secondary reaction 
channels comes from the use only of the uncorrelated part of the nucleon spectral function. Therefore, inclusive 
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subthreshold kaon production in pA collisions does not provide an information on the high momentum components 
within target nucleus. 
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